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The anharmonicity of a magnetic atomic trap at long distances from its center (forty times the 
radius of the first atomic orbit along the elongated axis, i.e. hundreds of microns) is shown to 
increase up to 45 % the temperature for Bose-Einstein condensation. This effect is perhaps small 
in the traditional traps, but should certainly be taken into account in magnetic microtraps, which 
characteristic dimensions are of order 1 mm. 

PACS numbers: 03.75.Fi 



Unlike the situation in liquid Helium, atomic vapors 
undergoing Bose-Einstein condensation (BEG) in mag- 
netic traps are very rarefiedi. As a consequence, the 
free-boson model works extremely well. Interaction ef- 
fects, computed in mean-field approximation, are shown 
to decrease the critical temperature for BEG in a few 
percents^-. Still lower corrections come from the approx- 
imation of the one-particle discrete spectrum by a con- 
tinuum of states, the so-called "finite-TV" corrections^, 
which also decrease the critical temperature. 

In the present paper, we show that a relatively im- 
portant increase of the critical temperature (up to 45 
%) could be related to trap anharmonicities at "long" 
distances from the center. By long distances, we mean 
around forty times the radius of the first atomic orbit 
along the elongated direction in the trap potential, i.e. 
around 200 /im. In the commonly used traps, which 
characteristic dimensions are a few centimeters, the be- 
lief is that anharmonic effects should be very weak at 
distances of 200 /xm. Nevertheless, we present an exam- 
ple of a recent experiment- , in which the measured values 
of the total number of atoms in the trap and the critical 
temperature are not consistent unless anharmonicity (or 
other) effects raising Tc are included. On the other hand, 
the trap potential should certainly be anharmonic at 200 
/im in the recently developed atomic microtraps^, which 
characteristic dimensions are of order 1 mm. 

In the continuum free-boson model, the critical tem- 
perature is simply estimated from the conditions A'o — 0, 
/i = in: 
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iVo is the number of atoms in the condensate, and e > is 
the one-atom excitation energy. The chemical potential 
takes values /i < 0. g is the density of states. 

Magnetic traps with cylindrical symmetry are very 
common. The trap potential, coming from the interac- 
tion of an hyperfine atomic species with the magnetic 
field, is proportional to the latter which, near the trap 
center, has a minimum and is written as: 



B{p,z) = Bo + bpp^ + b, 
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The corresponding harmonic oscillator frequencies are cup 
and LUz- It is useful to define cvq = (w^Wz)^/"^, in terms 



of which the level density is written g(e) 
and the critical temperature becomes: 
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where C(a;) = X]^5°=i 1/?^^ is the Riemann Zeta function. 

The corrections to Tc due to interaction or finite-A^ 
effects take, respectively, the formSi^: 
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where ro = ^/K/Jmuooj, a is the s-wave scattering length, 
and uj = {ujz + 2ujp)/3. Typical values N ^ 10^, 
a/ro ^ 10~^ lead to small numbers in the r.h.s. of Eqs. 
(|4I5() . As mentioned, the corrections I0J are computed 
within mean-field theory. More elaborated variational 
calculations based on diffusion Monte-Garlo^ or hyper- 
netted chain theory^ support the validity of mean-field 
approximations at the relevant particle densities. Notice 
that both corrections Q and Q are small and make Tc 
decrease. In a general power-law trap, the corrections 
may be of any sign and their relative magnitudes may 
be still greater^. However, for usual magnetic traps one 
expects a harmonic potential coming from at least 
in the vicinity of the trap center. 

Let us suppose that at long distances, r ^ rg, anhar- 
monic terms are relevant in In fact, one expects a 
saturation of the magnitude of B, which should approach 
the bias field value as the distance increases. In a trap 
potential that saturates at long distances, the density of 
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FIG. 1: Shift in the critical temperature as a function of the 
energy cutoff. Solid line: corrections coming from Eq^ 
crosses: finite- calculations for the trap used in Ref. 



energy levels decreases at high energiesS-. It means that 
the occupation of the ground state is increased, and thus 
the critical temperature is raised. 

We will model the density of levels of the (saturating) 
anharmonic potential in the simplest way: by truncating 
the one-particle spectrum. It means that an upper inte- 
gration limit, emax, is set in (Q. The critical temperature 
is thus determined from: 



t\{3) = C(3)-ii3(e-*")-teLi2(e~*") 
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where we have defined t = /Tc, x = emax/{kBT^), and 
Lik{z) — z"^ /n^ is the polylogarithm function. 

The corresponding correction to T^^, i. e. 1/i — 1, is 
shown in Fig. ^ as a function of x. It is around 10% 
when €max ~ 'iksT^- As ksT^ ~ IQ^fiWz in many of 
the experiments, energies of order AksT^ involve particle 
orbits of radius around 60 tq^ . On the other hand, when 
^max/kB ~ 2.5 T°, corrections to T° reach 40 %, and the 
cutoff distance is around 40 roz. A fit to the long-x tail 
of the curve leads to: 
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In quality of example, let us consider the trap used 
m Ref. 13, in which BEG of spin-polarized Helium was 
achieved. The example is quite interesting because inde- 
pendent measurements of N and Tc are reported. The 
authors-provided frequencies are i^z — 115 Hz, i^p — 1090 
Hz. The critical temperature, number of atoms in the 
trap and the scattering length are estimated as Tc = 
4.7 ± 0.5 fiK, N = {5± 2.5) x 10^ a « 16 nm, respec- 
tively. Turning back to Eqs. (|3I5|I . we obtain that a 
temperature Tc — 4.7 fiK, corresponds to a number of 
atoms N — 1.4 x 10^, well above the error bars. On 
the other hand, if we take the lowest value Tc — 4.2 fiK, 
then the number of atoms is N = 8 x 10^, a value closer 
but still outside error bars. The reason for this apparent 
discrepancy may be either the rough estimation of TV, 
or the anharmonicity corrections, which have not been 
included. 

In conclusion, we have shown that the saturation of the 
magnetic field at distances around 200 - 300 fim from the 
trap center leads to a decrease of the density of energy 
levels, and thus to an increase of the ground-state oc- 
cupation and to an increase of Tc- We have modeled 
the decrease of the level density simply by means of an 
energy cutoff in the single-particle spectrum. The cor- 
responding correction to T^ partially cancels out with 
interaction or finite-iV effects, both of which lower the 
critical temperature. Although expected to be small, de- 
tailed calculations of the magnetic field configurations in 
the commonly used traps are needed in order to evaluate 
the role of anharmonicity corrections. With regard to 
the atomic microtraps, which characteristic dimensions 
are of order 1 mm, anharmonicities are expected to be 
very important and should be taken into account. 
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